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Abstract — Micromachining  techniques  are  being  increasingly 
used  to  develop  miniaturized  sensor  and  actuator  systems.  These 
system  designs  tend  to  be  captured  as  layout,  requiring  extraction 
of  the  equivalent  microelectromechanical  circuit  as  a  necessary 
step  for  design  verification.  This  paper  presents  an  extraction 
methodology  to  (re-)construct  a  circuit  schematic  representation 
from  the  layout,  enabling  the  designer  to  use  microelectrome¬ 
chanical  circuit  simulators  to  verify  the  functional  behavior  of  the 
layout.  This  methodology  uses  a  canonical  representation  of  the 
given  layout  on  which  feature-based  and  graph-based  recognition 
algorithms  are  applied  to  generate  the  equivalent  extracted 
schematic.  Extraction  can  be  performed  to  either  the  atomic 
level  or  the  functional  level  representation  of  the  reconstructed 
circuit.  The  choice  of  level  in  hierarchy  is  governed  by  the  trade  off 
between  simulation  time  and  simulation  accuracy  of  the  extracted 
circuit.  The  combination  of  the  MEMS  layout  extraction  and 
lumped-parameter  circuit  simulation  provides  MEMS  designers 
with  VLSI-like  tools  enabling  faster  design  cycles,  and  improved 
design  productivity.  [682] 

Index  Terms — Atomic  elements,  canonical  representation, 
extraction,  functional  elements,  microelectromechanical  systems 
(MEMS). 


I.  Introduction 

ICROELECTROMECHANICAL  systems  (MEMS) 
integrating  multidomain  sensors  and  actuators  using 
conventional  batch  microfabrication  processes  are  becoming 
increasingly  complex.  Manual  verification  of  layouts  of 
such  complex  systems  is  virtually  impossible.  This  has  led 
to  an  increasing  need  for  MEMS  layout  verification  tools. 
Conventional  verification  tools  using  finite  element  (FEM) 
analysis  or  boundary  element  (BEM)  analysis  tend  to  be  quite 
cumbersome  and  time  consuming  for  such  large  designs.  This 
paper  proposes  an  alternative  approach  of  using  an  extractor 
[1]  which  reads  in  the  geometric  description  of  the  layout  and 
reconstructs  the  corresponding  MEMS  circuit.  This  enables  the 
designer  to  use  the  circuit-level  lumped  parameter  simulators 
[3] — [6]  for  faster  and  more  convenient  layout  verification.  In 
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addition,  minor  errors  in  the  layout,  like  missing  connections, 
can  be  easily  detected  in  the  reconstructed  circuit. 

Extraction  is  common  in  the  VLSI  world  where  the  main 
challenge  lies  in  extracting  parasitic  capacitances  and  resis¬ 
tances  of  long  interconnect  wires  [7] — [9].  In  the  MEMS  world 
we  have  similar  challenges  along  with  added  consideration  to 
geometrical  features  like  orientation,  relative  location,  etc., 
from  a  mechanical  perspective.  For  example,  a  L-shaped 
structure  may  be  recognized  as  an  interconnect  with  lumped 
or  distributed  parasitic  capacitance  and  resistance  by  a  VLSI 
extractor;  but,  depending  on  the  context,  a  MEMS  extractor 
might  recognize  it  as  a  crab  leg  spring.  Thus,  MEMS  extraction 
needs  to  inherit  the  principles  of  VLSI  extraction  [10],  [11]  and 
also  expand  on  them  by  incorporating  geometrical  heuristics 
from  the  mechanical  world.  Traditional  VLSI  extraction  tools 
are  designed  to  recognize  only  electrical  elements  like  transis¬ 
tors,  resistances,  capacitances.  They  do  not  currently  support 
element  recognition  for  mechanical  elements  like  beams, 
plates,  joints  which  are  essential  in  MEMS.  However,  the  basic 
geometry  processing  functions  needed  for  the  extraction  of 
mechanical  and  electrical  elements,  like  proximity  analysis, 
shape  analysis,  area  calculation,  etc.,  are  the  same.  Access  to 
such  geometry  processing  functions  would  help  us  modify 
existing  VLSI  extraction  tools  to  include  the  element  recog¬ 
nition  modules  needed  for  extracting  mechanical  elements  in 
addition  to  the  electrical  elements.  Since  none  of  the  existing 
commercial  VLSI  extraction  tools  are  transparent  enough  to 
allow  access  to  their  core  geometry  processing  routines,  custom 
tools  are  necessary  for  MEMS  extraction.  This  paper  explains 
the  element  recognition  and  geometry  processing  algorithms 
written  to  implement  a  custom  prototype  extractor  capable  of 
extracting  MEMS  elements. 

In  VLSI,  designers  start  by  capturing  the  connectivity  and 
individual  element  parameters  in  the  form  of  a  schematic. 
Layout  in  VLSI  designs  merely  adds  parasitic  resistances 
and  capacitances  due  to  routing  and  placement  of  the  circuit 
elements.  These  effects  are  captured  in  the  extracted  view  of 
the  circuit  which  is  then  used  to  do  a  Layout  versus  Schematic 
(LVS)  check  [12]  for  layout  verification.  Unlike  the  VLSI 
world,  layout  forms  the  most  important  representation  for 
MEMS  designs  [13].  Here  the  placement  of  the  elements  is 
as  crucial  as  the  connectivity  and  parameters.  For  example,  a 
symmetrical  design  in  VLSI  world  would  simply  mean  a  sym¬ 
metrical  netlist  from  power  supply  to  ground,  while  in  MEMS  a 
symmetrical  design  would  additionally  mean  that  the  structures 
(like  springs  around  a  plate)  are  symmetric  geometrically.  This 
also  increases  the  importance  of  having  a  LVS  tool  for  MEMS. 
Initial  attempts  to  achieve  such  capability  is  limited  to  higher 
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level  connectivity  analysis  as  in  [14].  Here,  the  design  entry  is 
at  the  schematic  level  and  the  layout  is  generated  automatically. 
Each  of  the  elements  in  the  layout  are  tagged  with  the  corre¬ 
sponding  schematic  element  and  the  LVS  tool  verifies  the  pin 
to  pin  connectivity  of  this  tagged  layout.  Such  a  methodology 
does  not  allow  manual  layout  generation  and  also  fails  to 
capture  the  placement  related  errors  in  the  layout.  In  contrast, 
the  extractor  presented  here  can  be  used  as  an  initial  framework 
for  a  MEMS  LV S  tool  having  the  flexibility  of  both  schematic 
and  layout  design  entry.  The  extractor  results  in  a  netlist  which 
captures  the  connectivity  and  the  placement  of  the  actual  layout 
along  with  the  various  properties  of  the  recognized  elements. 
Such  a  netlist  can  be  easily  used  to  perform  a  LVS  check  as 
in  VLSI.  However  improvements  need  to  be  made  to  create 
a  complete  LVS  tool  capable  of  addressing  MEMS  specific 
errors  like  placement  errors,  geometrical  symmetry  errors, 
etc.  As  with  VLSI,  schematic  design  capture  and  automated 
layout  generators  will  eventually  replace  manual  layouts  for 
MEMS .  At  that  time,  extraction  will  still  be  needed  to  capture 
the  mechanical  parasitics  in  the  final  layout.  Examples  of 
mechanical  parasitics  include  joints  between  two  beams  of 
different  aspect  ratios  [15], 

This  paper  focuses  on  single  layer  suspended  MEMS 
components  because  of  their  wide  applicability.  This  class 
of  MEMS  components  cover  a  wide  range  of  applications, 
such  as  pressure  sensors  [16],  micromirrors  [17],  RF  switches 
[18],  accelerometers  [19]  and  resonator  filters  [20];  and  can 
be  fabricated  using  MEMS  fabrication  processes  like  bulk 
micromachining  [21],  LIGA  [22],  surface  micromachining 
[23]  or  DRIE  [21].  As  a  representative  process  we  will  use  the 
Multi-User  MEMS  Process  (MUMPS)  [24]  which  is  similar 
to  many  other  surface  micromachining  processes  like  Sandia’s 
SUMMiT  process  [25],  Analog  Devices’s  iMEMS  process  [26] 
and  Case  Western  Reserve  University’s  Polycrystalline  Silicon 
Carbide  surface  micromachining  process  [27],  An  example  of 
a  MEMS  component  built  using  this  process  is  shown  in  Fig.  1. 
It  consists  of  an  I-shaped  plate  suspended  by  a  pair  of  folded 
flexure  springs.  Two  sets  of  electrostatic  comb  drives  on  either 
side  act  as  electromechanical  transducers.  The  process  uses 
two  layers  of  polycrystaline  silicon  separated  by  a  sacrificial 
oxide  layer.  The  oxide  is  etched  away  in  the  final  step  of  the 
process  releasing  the  suspended  mechanical  structures  (defined 
by  the  structural  polysilicon  mask )  which  are  connected  to  the 
base  conducting  polysilicon  layer  at  the  anchors  (defined  by 
the  anchor-cut  mask).  Additional  holes  (defined  by  the  hole 
mask )  are  used  on  the  large  areas  of  the  structural  layer  to  aid 
the  release  step.  These  areas  also  require  bushings  (defined 
by  the  dimple  mask)  to  prevent  sticking  to  the  bottom  surface 
during  the  release  step.  Though  the  current  extraction  tool  uses 
the  mask  conventions  of  MUMPS  for  the  recognition  process, 
it  can  be  very  easily  extended  to  other  processes. 

A  MEMS  component  can  be  hierarchically  decomposed  into 
functional  elements  like  springs  and  comb  drives.  Each  of  these 
functional  elements  can  in  turn  be  decomposed  into  more  funda¬ 
mental  or  atomic  elements  like  beams,  joints,  fingers,  anchors, 
plates  and  gaps.  The  various  atomic  and  functional  elements  for 
the  example  in  Fig.  1  are  marked  on  the  figure.  A  MEMS  com¬ 
ponent  can  be  represented  using  schematic  symbols  and  asso- 


mmmnmnn 

ouitmtuHitui 


anchor-cut 


hushing 
W Y  (dimple) 
\  structural 
polysilicon 
interconnect 
polysilicon 


plate 


folded 

flexure 

spring 


anchor 

finger 

gap 


comb 

drive 


beam 


(a) 


(b) 


Fig.  1 .  (a)  A  scanning  electron  micrograph  (SEM)  picture  of  a  folded  flexure 

comb  drive  micro  resonator  fabricated  in  the  MUMPS  process  with  atomic 
elements  (listed  on  the  right  of  the  SEM)  and  functional  elements  (listed  on  the 
left  of  the  SEM),  (b)  3-D  view  of  the  area  highlighted  in  (a). 


dated  behavioral  models  at  each  hierarchical  level  resulting  in 
either  an  atomic  level  or  a  functional  level  circuit  schematic. 
The  functional  level  schematic  simulates  faster  than  the  atomic 
level  schematic.  The  extraction  tool  gives  the  user  the  freedom 
to  choose  the  hierarchy  level  depending  on  the  specific  needs  of 
accuracy  and  simulation  time. 

The  sections  that  follow  give  a  detailed  description  of  the 
algorithms  and  representations  used  in  the  prototype  extractor. 
The  extraction  process  is  divided  into  three  steps.  First,  an 
unique  representation  for  the  given  layout  geometry,  referred 
to  as  the  canonical  representation,  is  generated.  Section  II 
motivates  the  need  for  such  a  representation  and  also  describes 
the  algorithm  to  create  a  canonical  representation.  The  next 
step  is  extraction  of  atomic  elements  which  is  described  in 
Section  III.  The  user  may  generate  an  atomic  level  schematic 
at  this  stage  or  proceed  toward  extraction  to  the  functional 
level.  Section  IV  describes  algorithms  used  to  detect  functional 
elements  like  springs  and  comb  drives.  Section  V  presents 
results  demonstrating  the  usefulness  of  the  extractor  followed 
by  the  conclusion  in  Section  VI. 

II.  Canonical  Representation 

Designers  tend  to  design  layouts  in  their  own  specific  way 
resulting  in  nonunique  representations  of  the  same  layout.  The 
extraction  tool  converts  such  representations  to  an  unique  rep¬ 
resentation  in  order  to  simplify  the  recognition  algorithms  used 
later  for  extraction.  There  are  various  types  of  unique  represen¬ 
tations  that  are  followed  in  VLSI  CAD,  out  of  which  the  tile 
plane  representation  [28]  is  well  known.  The  canonical  repre¬ 
sentation  is  a  derivative  of  the  tile  plane  representation.  Un¬ 
like  the  tile  plane  representation,  where  each  tile  can  have  mul¬ 
tiple  neighbors  on  each  of  its  sides,  the  polygons  in  canon¬ 
ical  representation  can  either  have  one  or  no  neighbors  on  each 
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edge.  Having  such  an  unique  neighbor  information  makes  the 
neighbor  based  recognition  algorithms  simpler  to  implement. 

As  a  vast  majority  of  MEMS  layouts  are  Manhattan,  the  pro¬ 
totype  implementation  is  limited  to  Manhattan  designs  only. 
Fig.  2  explains  canonical  representation  in  the  context  of  Man¬ 
hattan  designs.  We  define  the  canonical  representation  of  the 
layout  to  be  the  one  which  uses  minimum  number  of  rectangles 
to  cover  the  given  layout  area,  such  that  infinitesimal  outward 
extensions  of  an  edge  of  any  rectangle  never  intersects  with  the 
interior  of  the  layout  area.  We  use  the  term  layout  area  to  define 
the  area  which  represents  the  actual  component  in  the  layout, 
i.e.,  it  is  the  interior  area(s)  defined  by  the  boundary  /boundaries 
of  the  geometrical  representation  of  the  component  in  the  layout. 
Thus,  in  the  canonical  representation,  the  Manhattan  layout  is 
made  up  of  rectangles  such  that  each  rectangle  has  at  most  one 
neighbor  per  edge  and  each  edge  is  either  fully  covered  by  a 
neighbor  or  not  covered  at  all.  This  can  be  easily  achieved  by 
extending  the  boundary  edges  into  the  interior  of  the  layout  area 
till  it  meets  another  boundary  edge.  The  resulting  representation 
uniquely  partitions  the  layout  area.  Extension  to  non-Manhattan 
geometries  can  be  accomplished  by  following  the  same  prin¬ 
ciple  after  replacing  rectangles  by  polygons. 

The  process  of  canonization  starts  from  the  input  geometric 
description  of  the  chip  (written,  perhaps,  in  CIF,  i.e.,  Caltech 
Intermediate  Form).  The  hierarchical  layout  description  of  the 
chip  is  first  flattened  followed  by  the  removal  of  all  overlaps 
between  polygons.  An  initial  rectangular  cover  of  the  layout  is 
then  obtained  for  the  suspended  structural  pattern  in  the  struc¬ 
tural  polysilicon  mask.  This  serves  as  input  to  the  actual  canon¬ 
ization  routine. 

The  primary  interaction  in  the  canonization  process  takes 
place  between  two  sets;  the  input  set  and  the  output  set.  The 
output  set  is  always  kept  in  canonical  state  with  respect  to  its 
contents  and  will  eventually  contain  the  canonical  version  of  the 
input  set.  Elements  from  the  input  set  are  selected  sequentially 
and  added  to  the  output  set.  Whenever  there  is  an  addition  to  the 
output  set,  its  equilibrium  might  be  destroyed  (i.e.,  the  output 
set  might  no  longer  be  a  canonical  set).  If  this  occurs,  a  series 
of  operations  is  initiated  which  ultimately  brings  the  output  set 
back  to  its  equilibrium  or  canonical  state.  This  is  repeated  till 
the  input  set  is  empty  at  which  point  the  output  set  will  contain 
the  canonical  representation  of  the  input  layout.  The  process 
which  drives  the  output  set  to  equilibrium,  after  it  is  disturbed 
by  the  insertion  of  a  new  element,  is  described  below: 

ADJACENT (rectangle  a,  rectangle  b)  :  re¬ 
turns  TRUE  if  a  and  b  have  a  common  edge 
SPLIT ( rectangle  a,  rectangle  b)  :  splits 
a  by  edges  of  b,  if  any  of  the  vertices 
of  the  edges  of  b  lies  on  any  edge  of  a, 
and  also  updates  the  neighbor  informa¬ 
tion 

CANONIZE (rectangle_set  R) 

G  =  NULL 
while  R\  —  NULL 

r  —  pop[R] 

P  =  NULL 

push[P,  r] 


Various  representations  of  the  same  structure 


tn  =  top  neighbor 
bn  =  bottom  neighbor 
rn  =  right  neighbor 
In  =  left  neighbor 
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Fig.  2.  Fully  partitioned  canonical  representation. 


Q  —  {a:  [ADJACENT (a:,  r);  x  is  an  element  of 
G} 

G  =  G  -  Q 

for  i  —  0  to  length[Q] 
for  j  =  0  to  length[P] 

SPLIT (P [ j ] ,  Q [ i ] ) 
for  i  —  0  to  length[P] 
for  j  —  0  to  length[Q\ 

SPLIT  (Q\j\,  P[i]) 

M  -  Q 

until  M  =  NULL 
N  =  NULL 

for  i  —  0  to  length[M] 

S  —  {a:|a:  is  neighbor  of  M[i];  x  is 

element  of  G} 

G  =  G  -  S 

for  j  —  0  to  length[S] 

SPLIT  (A  ./j.  M  i\) 

N  =  N  +  S 
P  =  P  +  M 
M  =  N 
G  =  G  +  P 

return 

The  neighbor  information  of  each  rectangle  in  the  resulting 
canonical  representation  comes  as  a  by-product  of  the  algo¬ 
rithm.  Each  rectangle  in  this  representation  has  four  pointers 
which  point  to  the  neighbor,  if  it  exists,  on  each  of  the  sides 
(as  shown  in  Fig.  2).  The  algorithm  has  a  worst  case  asymp¬ 
totic  upper  bound  of  0(n1-0)  (See  the  “Appendix”),  where  n  is 
the  number  of  rectangles  in  the  final  canonized  representation. 
Though  this  is  much  greater  than  the  usual  complexity  of  VLSI 
CAD  tools,  which  normally  tend  to  be  0(nlgn ),  it  does  not  pose 
a  serious  time  restriction  since  the  problem  size  is  much  smaller 
(~1000)  than  that  encountered  in  VLSI  world.  Nevertheless,  the 
algorithm  does  have  ample  scope  of  improvement  in  terms  of 
time  and  storage  requirements  if  the  problem  size  increases. 

III.  Extraction  of  Atomic  Elements 

Atomic  elements  [as  listed  on  the  right  side  of  Fig.  1(a)]  form 
the  fundamental  building  blocks  of  MEMS  components  [29]. 
This  section  describes  some  common  atomic  elements  that  are 
used  to  build  suspended  MEMS  components  followed  by  the 
description  of  their  recognition  heuristics  and  algorithms. 
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A.  Atomic  Elements 

Suspended  components  essentially  consists  of  structural 
areas  suspended  over  the  substrate.  Such  suspended  areas  can 
be  partitioned  into  two  groups  based  on  their  relative  rigidity: 
plates  and  beams.  The  suspended  areas  of  the  structure,  which 
are  relatively  rigid  to  forces  along  the  plane  of  the  structure, 
are  defined  to  be  plates.  They  are  the  major  contributors  to  the 
mass  of  the  suspended  structure.  The  structural  compliance  of 
the  suspended  structure  is  decided  by  the  beams.  Geometrically 
these  are  rectangular  areas  having  neighbors  only  on  their 
shorter  sides.  Cantilever  beams  are  often  classified  separately 
as  fingers  and  are  extensively  used  to  design  electrostatic 
actuators  and  sensors.  Sometimes  fingers  are  provided  with 
pedestals  to  reduce  the  inter-finger  gap  below  the  lithography 
limit  for  higher  sensitivity  of  the  mechanical  to  electrical 
transfer  function  for  the  electrostatic  transducer  [30],  Two  or 
more  beams  are  connected  by  joints  which  can  be  modeled 
in  the  adjacent  beams  for  simplicity.  Hence  joints  serve  as 
logical  connectivity  elements.  The  suspended  structure  is 
connected  to  the  base  (interconnect  polysilicon)  at  the  anchors 
which  are  defined  by  the  anchor-cut  mask.  These  areas  provide 
electrical  connection  to  the  suspended  structure  and  also  act  as 
mechanical  pillars  supporting  the  suspended  areas. 

B.  Extraction  Flow  for  Atomic  Elements 

At  this  stage  of  extraction,  a  canonical  representation  of  the 
layout  geometry  is  available  and  the  objective  is  to  tag  each 
rectangle  in  the  canonical  representation  by  the  correct  element 
type  (anchor,  plate,  finger,  beam  or  joint).  The  recognition  of 
each  type  of  elements  is  a  two  step  process  as  shown  in  the  flow 
diagram  in  Fig.  3.  The  first  step  marks  probable  elements  which 
are  confirmed  in  the  second  step  using  stricter  rules. 

The  first  step  in  the  recognition  process  is  to  mark  potential 
plate  and  anchor  areas  using  information  from  the  nonstructural 
masks.  For  example,  overlap  of  the  structural  mask  with  the  an¬ 
chor-cut  mask  is  used  to  tag  potential  anchors.  Similarly  overlap 
with  hole  mask  and  bushing  mask  are  used  to  tag  potential  plate 
rectangles. 

The  next  step  is  to  mark  potential  fingers  from  the  neighbor 
information  of  the  rectangles.  These  are  rectangles  having  a 
neighbor  at  one  of  its  shorter  sides.  The  subroutine  to  confirm 
fingers  selects  those  rectangles  or  connected  sets  of  rectangles 
(fingers  split  into  multiple  rectangles  during  the  canonization 
process),  out  of  the  set  of  potential  fingers,  which  satisfy  the 
criteria  that  a  finger  can  have  only  one  neighbor  and  it  should 
lie  only  on  one  of  its  shorter  side.  If  fingers  are  detected  they 
are  removed  from  the  total  set  of  rectangles  and  the  remaining 
structure  is  recanonized.  This  reduces  the  number  of  rectangles 
being  checked  by  the  subroutines  and  also  helps  in  detecting 
pedestals  which  come  out  as  fingers  in  the  remaining  geometry. 
The  worst  case  complexity  of  the  finger  detection  is  0{n  +  mp) 
where,  n  is  the  initial  number  of  rectangles,  to  number  of  rect¬ 
angles  tagged  in  the  first  loop  and  p  is  the  remaining  rectangles 
(n  —  to).  Though  theoretically  this  implies  quadratic  time  com¬ 
plexity,  in  reality  it  does  not  create  any  serious  time  limitation 
because  the  problem  size  is  small  (to  ~  100). 

In  the  next  step,  any  rectangular  empty  space  surrounded  by 
filled  rectangles  on  each  side  is  tagged  as  a  potential  hole.  This 
is  achieved  by  sieving  out  rectangles  in  the  canonized  represen- 


Fig.  3.  Flow  chart  for  extraction  of  atomic  elements. 

tation  of  the  geometric  not  of  the  structural  geometry.  This  is 
followed  by  the  detection  of  potential  beams  which  are  rectan¬ 
gles  having  neighbors  on  each  of  their  shorter  sides.  This  is  used 
by  the  hole  confirmation  routine  which  removes  potential  holes 
that  are  actually  gaps  between  beams.  The  holes  detected  are 
then  replaced  by  filled  rectangles  followed  by  a  recanonization 
of  the  resulting  geometry.  The  change  in  physical  parameters 
(like  mass-factor,  center  of  mass,  etc.)  due  to  the  addition  of 
these  filled  rectangles  is  annotated  to  the  corresponding  rect¬ 
angles.  In  addition  to  reducing  the  number  of  rectangles  in  the 
representation  of  the  layout,  the  recanonization  also  reduces  the 
chances  of  split  beams. 

The  next  step  detects  beams  and  joints  using  neighbor-based 
heuristics  followed  by  the  final  detection  of  plates  and  anchors. 
Previously  marked  potential  plate  and  anchors  are  used  as  seeds 
to  recursively  expand  into  unrecognized  areas  and  mark  them  as 
plates  or  anchors.  The  expansion  process  checks  each  rectangle 
only  once  and  hence  runs  in  linear  time.  The  final  recognized  set 
of  rectangles  can  be  used  to  generate  an  atomic  level  schematic 
or  can  be  used  to  proceed  to  functional  level  extraction. 

IV.  Extraction  of  Functional  Elements 

Complex  MEMS  components  are  best  modeled  in  the  func¬ 
tional  level  schematic.  Such  schematics  have  fewer  elements 
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than  atomic  level  schematics  and  hence  allow  faster  simulations 
without  any  significant  change  in  accuracy  of  result  [31].  This 
motivates  the  extraction  to  the  functional  level  instead  of  stop¬ 
ping  at  the  atomic  level.  This  section  describes  the  various  types 
of  functional  elements  [as  listed  on  left  side  of  Fig.  1(a)]  and  the 
algorithms  that  are  used  for  extracting  them. 

A.  Plate  and  Anchor 

The  canonization  process  results  in  a  large  number  of  in¬ 
terconnected  plates  and  anchors  in  the  atomic  level  schematic. 
Combining  such  interconnected  rigid  plates  and  anchors  to  get 
a  minimal  representation  would  result  in  a  significant  decrease 
in  the  number  of  elements  in  the  schematic.  The  plates  and  an¬ 
chors  in  such  a  minimal  representation  are  referred  to  as  func¬ 
tional  plates  and  anchors.  The  algorithm  described  below  is  used 
to  achieve  a  maximal  horizontal  or  maximal  vertical  representa¬ 
tion  of  interconnected  sets  of  similar  atomic  elements  depending 
on  which  representation  results  in  smaller  number  of  elements. 

HORIZONTAL ( rectangle_set  A):  returns  the 
maximal  horizontal  representation  of  A 
VERTICAL (rectangle_set  A)  :  returns  the 
maximal  vertical  representation  of  A 
PATH (rectangle  a,  rectangle  b)  :  returns 
TRUE  if  a  path  exists  from  a  to  b  such 
that  every  element  of  the  path  including 
a  and  b  are  of  the  same  type 
OPTIMIZE_ELEM (rectangle_set  G ) 

H  =  NULL 
V  =  NULL 
while  G\  =  NULL 

9  =  Pop[G\ 

P  =  NULL 
push[P,  <7] 

Q  —  {:r|PATH(<7,  x);  x  is  an  element  of  G} 

P  -  P  +  Q 
G  —  G  -  Q 

A!  =  HORIZONTAL  (P) 

A  =  VERTICAL  (A) 

B'  =  VERTICAL  (P) 

B  =  HORZONTAL(H') 

H  —  H  +  A 
V  =  V  +  B 

if  length[H]  <  length[V]  then  return(if) 
else  return(V) 

Since  the  neighbor  information  for  each  rectangle  is  already 
available  from  the  canonization  algorithm,  the  sets  Q  can 
be  easily  obtained  by  visiting  each  rectangle  only  once.  The 
merging  algorithms  also  make  use  of  the  neighbor  information 
and  hence  check  each  rectangle  in  the  set  only  once.  Hence  the 
entire  algorithm  runs  in  linear  time. 

B.  Electromechanical  Comb  Actuators 

Silicon  microstructures  have  long  been  actuated  and  sensed 
electrostatically.  A  widely  used  electrostatic  actuator  is  the 
linear  comb  drive  [32]  made  up  of  interdigited  fingers  which 


/C 


s 


(c) 


Fig.  4.  Different  types  of  finger  arrangement  (a)  pair  of  cantilever  beams 
forming  the  building  block  of  electrostatic  comb  drive,  (b)  fingers  with 
pedestals,  (c)  differential  comb  finger  arrangement. 


Fig.  5.  Comb  drive  extraction  flow  chart. 


may  or  may  not  have  pedestals  [30]  [see  Fig.  4(a)  and  (b)].  One 
side  of  the  comb  is  fixed  (stator)  while  the  other  side  is  allowed 
to  move  (rotor).  Another  popular  electrostatic  actuator  is  the 
differential  comb  drive  [33]  [see  Fig.  4(c)]  using  three  sets  of 
electrically  isolated  comb  fingers  arranged  such  that  the  rotor 
set  of  comb  finger  sees  different  sets  of  capacitances  on  its  two 
sides.  Such  a  structure  is  used  to  sense  transverse  motion  via 
differential  sensing  of  the  two  sets  of  capacitances  and  has  an 
added  advantage  of  reduced  levitation  problems  [34]. 

The  comb  drive  extraction  flow  is  shown  in  Fig.  5.  It  starts 
with  a  connectivity  analysis  of  the  set  of  recognized  fingers. 
Electrically  connected  fingers  are  given  the  same  connectivity 
number  and  then  sorted  into  buckets  based  on  their  orientation 
and  connectivity.  Each  such  finger  bucket  is  then  checked  for 
uniformity  of  the  fingers  with  respect  to  their  geometrical  pa¬ 
rameters  like,  region  of  occurrence,  length  of  fingers,  width  of 
fingers  and  interfinger  gap.  If  the  fingers  have  pedestals,  then 
the  geometrical  parameters  of  the  pedestal  (like  region  of  oc¬ 
currence,  length  and  width  of  the  pedestal,  interpedestal  gap 
and  the  relative  position  of  the  pedestal  with  respect  to  the  thin 
cantilever  finger)  are  also  checked.  The  buckets  are  partitioned 
whenever  any  nonuniformity  is  found  in  any  of  these  param¬ 
eters.  A  set  of  overlapping  pair  or  triad  of  such  buckets  with 
different  electrical  connectivity  numbers  will  result  in  a  linear 
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Fig.  6.  Springs  in  the  spring  library:  (a)  fixed-fixed,  (b)  crab  leg,  (c)  U-spring, 
(d)  meander  spring,  and  (e)  folded  flexure. 


or  differential  comb  drive  respectively.  The  box  covers  of  the 
buckets  are  checked  using  box  overlap  rules  to  find  such  pairs 
or  triads.  Any  such  set  is  then  checked  to  avoid  uncoupled  fin¬ 
gers  and  finally  grouped  together  to  form  a  linear  or  differential 
comb  drive. 

C.  Mechanical  Springs 

Springs  are  composed  of  beams  and  joints  and  connect  the 
suspended  plate  to  the  anchors.  Few  commonly  used  springs 
are  shown  in  Fig.  6.  The  fixed-fixed  flexure  [see  Fig.  6(a)]  con¬ 
sists  of  a  simple  straight  beam  connecting  the  suspended  plate 
to  the  anchor  and  has  a  high  spring  constant  because  of  exten- 
sional  axial  stress  in  the  beams.  Crab-leg  springs  and  U-springs 
[see  Fig.  6(b)  and  (c)]  are  modifications  of  the  fixed-fixed  beam 
to  reduce  peak  stress  in  the  flexure  at  the  cost  of  reduced  stiff¬ 
ness  in  undesired  directions.  A  meander  spring  [see  Fig.  6(d)] 
is  also  a  modified  version  of  fixed-fixed  flexure  which  helps 
achieve  more  compliance  using  less  space.  A  folded  flexure  [see 
Fig.  6(e)]  design  reduces  axial  stress  and  also  has  the  advantages 
of  providing  more  compliance  while  occupying  less  area.  The 
springs  types  (like  crab  leg,  serpentine  spring,  folded  flexure, 
etc.)  to  be  recognized  are  stored  in  the  form  of  graphs  in  a  li¬ 
brary  file.  The  spring  detection  routine  reads  in  this  library  file 
to  create  an  internal  copy  of  the  graphs.  Any  spring  consisting 
of  a  contiguous  set  of  beams  and  joints  that  can  be  represented 
in  the  manner  described  below  can  be  recognized  by  the  spring 
detection  routine.  However,  the  methodology  described  can  be 
generalized  to  handle  springs  made  of  other  atomic  elements  by 
incorporating  user  definable  atomic  elements. 

The  spring  detection  routine  stores  the  spring  library  by  cre¬ 
ating  a  finite  state  machine  (FSM)  for  each  of  the  springs  de¬ 
fined  in  the  library.  Each  of  the  FSMs  can  be  defined  by  M  = 
{Q,  S,  L,  G,  F,  X},  where 

Q  :  states  =  [S,  [intermediate  states],  F,  X}; 

S  :  start  state  =  anchor  point; 

L  :  inputs  =  { [joints],  [beams],  NULL}; 

G  :  transition  rules; 


TABLE  I 

Dictionary  of  Joints 


Joint 

name 

m  - 

param 

t- 

param 

ports 

example 

j+ 

+i 

0 

2 

j. 

-i 

0 

2 

JT0 

0 

0 

3 

JT+ 

+1 

+  1 

3 

Jr. 

-1 

+  1 

3 

“V 

Jo 

0 

+1 

4 

4-* 

layout  with  atomic 
elements 


make  contiguous  sets 
of 

beams  and  joints 


spring 

detected 


) 


no  (  undetected 
spring 


D 


Fig.  7.  Spring  extraction  flow  chart. 

F  :  set  of  final  states;  and 

X  :  exit  state. 

A  joint,  in  such  a  language,  is  defined  to  be  a  node  having  one 
input  port  and  at  most  three  output  ports  and  is  labeled  using 
the  “to”  (from  moment)  and  “t”  (from  transition)  parameters. 
The  ^-parameter  is  1  only  if  there  is  an  output  port  along  the 
direction  of  the  input  port.  An  output  port  at  right  angles  to  the 
input  port  contributes  a  +1  or  —  1  to  m-parameter  depending  on 
the  direction  of  (anticlockwise  or  clockwise  respectively).  The 
six  types  of  joints  possible  using  such  a  convention  are  shown 
in  Table  I.  The  set  of  beams  for  the  language  depend  on  the 
spring  to  be  detected.  For  example,  a  crab  leg  spring  requires 
two  beams  [see  Fig.  6(b)]  which  may  or  may  not  be  equal  in 
dimension,  while  a  folded  flexure  requires  four  type  of  beams 
which  must  be  arranged  as  shown  in  Fig.  6(d). 

Connected  sets  of  beams  and  joints  obtained  after  the  atomic 
recognition  are  passed  through  each  of  these  FSM’s  to  recog¬ 
nize  their  type.  For  each  such  set,  the  input  is  started  from  a 
beam  which  is  connected  to  an  anchor  rectangle.  The  flow  of 
the  spring  detection  algorithm  is  shown  in  Fig.  7.  Any  unde¬ 
tected  spring  is  replaced  by  its  atomic  level  representation,  i.e., 
in  the  form  of  beams  and  joints. 
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beam(a) 


beam(a) 


beam(b) 

joint(J+) 

joint(J.) 

beam(b) 


M  =  {Q,  S,  L,  G,  F,  X) 
Q  =  {S,1,2,3, 4,  F,  X} 
L  =  {a,  b,  J+,  J_,  <}>} 

G  =  transition  rules 
F  =  final  state 
X  =  exit  state 


Fig.  8.  FSM  to  recognize  a  crab  leg. 


Fig.  8  shows  the  FSM  that  defines  a  crab  leg.  For  a  connected 
set  of  beams  and  joints,  the  algorithm  first  detects  the  beam  that 
is  connected  to  an  anchor  rectangle  and  initializes  the  current 
state  to  be  the  start  state  (S  in  Fig.  8)  of  the  FSM.  At  every 
state  the  algorithm  uses  the  next  element  in  the  connected  set 
of  beams  and  joints  as  the  input  to  decide  the  next  state.  For 
example,  for  the  start  state,  if  the  next  input  is  a  beam  (which, 
being  the  first  instance  of  beam,  will  be  registered  as  beam  a) 
then  the  current  state  advances  to  state  1.  If  the  input  is  not  a 
beam  (i.e.,  if  it  is  a  joint  or  if  there  are  no  other  elements  in  the 
set)  then  the  current  state  advances  to  state  x.  Reaching  state  x 
implies  that  the  given  set  is  not  a  crab  leg  and  the  algorithm  exits 
from  the  current  spring  type  (crab  leg).  It  then  tries  to  match 
the  set  with  some  other  spring  type  in  the  library.  If  there  are 
no  more  spring  types  remaining  in  the  library,  the  current  set  is 
marked  as  an  unrecognized  spring.  For  a  set  that  forms  a  crab 
leg,  the  first  input  will  be  a  beam  ( beam  a)  followed  by  either  a 
J_l_  or  J_  joint  and  finally  another  beam  (which  will  be  marked 
as  beam  b).  At  this  point  the  current  state  will  reach  state  4.  For  a 
crab  leg,  there  should  be  no  more  elements  (beams  or  joints)  left 
as  input.  If  this  is  the  case  (i.e.,  the  input  is  NULL,  symbolized 
by  in  Fig.  8),  then  the  current  state  reaches  final  state  (F)  and 
the  set  of  beams  and  joints  is  recognized  as  a  crab  leg. 

V.  Results 

This  section  presents  some  results  to  demonstrate  the  capa¬ 
bility  of  the  prototype  extractor  which  implements  the  algo¬ 
rithms  discussed  in  the  previous  sections. 

A.  Accelerometer 

Fig.  9  shows  the  results  for  an  accelerometer  design  [33] 
using  U  spring  as  the  flexure  and  a  pair  of  differential  comb 
drives  for  sensing  the  output.  Any  inertial  force  acting  on  the 
central  plate  due  to  an  external  acceleration  causes  the  sus¬ 
pended  structure  to  move  in  the  horizontal  direction  resulting 
in  a  differential  change  of  the  capacitances  in  the  comb  drive 
which  can  then  be  sensed  using  sensing  electronics.  Fig.  9(a) 
shows  the  initial  layout.  Notice  the  structural  holes  that  are 


(g) 


Fig.  9.  Accelerometer:  (a)  input  layout,  (b)  canonical  representation,  (c) 
intermediate  state,  (d)  atomic  elements  recognized,  (e)  functional  elements 
recognized,  (f)  reconstructed  functional  level  schematic,  and  (g)  transient 
response  to  a  1-g  pulse  acceleration  input. 


present  in  the  plate  area  of  the  layout.  Fig.  9(b)  shows  the 
fully  partitioned  canonical  representation  of  the  component. 
The  information  from  the  anchor-cut  mask  is  used  to  mark 
the  anchor  areas  of  the  layout  [see  Fig.  9(c)].  This  is  followed 
by  the  recognition  of  atomic  elements  [see  Fig.  9(d)].  The 
structural  holes  in  the  plate  have  been  replaced  by  actual 
plate  areas.  The  resulting  decrease  in  mass  is  annotated  in  the 
schematic  resulting  from  this  representation.  The  next  step  is  to 
recognize  the  functional  elements,  i.e.,  U  spring  and  differential 
comb  drive.  Fig.  9(e)  shows  the  final  recognized  state  of  the 
component.  The  corresponding  functional  level  schematic  is 
shown  in  Fig.  9(f).  The  transient  response  of  the  accelerometer 
due  to  a  pulse  acceleration  input  of  lg  is  shown  in  Fig.  9(g). 


BAIDYA  et  al.:  AN  EXTRACTION-BASED  VERIFICATION  METHODOLOGY  FOR  MEMS 


9 


error 

detected 

folded 

flexure 

spring 

simple 

comb 

drive 


(a)  (b)  can  be  simulated  to  verify  the  functional  behavior  of  the  device. 

This  example  demonstrates  the  capability  of  the  extractor  to  act 
Fig.  11.  Three-fold  symmetric  gyroscope:  (a)  input  layout  and  (b)  recognized  as  a  verification  t0ol  for  compiex  MEMS  devices. 


The  simulation  using  functional  level  schematic  was  found  to 
be  approximately  10  times  faster,  than  the  same  simulation 
using  atomic  level  schematic,  with  no  appreciable  change  in 
accuracy  of  the  result. 

B.  Erroneous  Resonator  Layout 

The  importance  of  the  extractor  is  demonstrated  in  this  ex¬ 
ample  (see  Fig.  10)  where  the  input  layout  of  a  folded-flexure 
resonator  [35]  was  found  to  have  a  very  small  error  which  was 
not  detected  by  the  human  eye.  When  extracted,  the  recognized 
representation  contained  two  sets  of  comb  actuators  instead  of 
just  one  pair.  On  inspecting  the  original  layout,  it  was  found 
that  there  was  a  difference  in  the  gaps  between  the  two  halves 
of  each  comb  drive.  This  was  because  when  the  half  was  being 
replicated  and  placed  to  double  the  size  of  the  comb  actuator, 
a  small  human  error  resulted  in  a  gap  which  was  more  than  the 
gaps  between  other  fingers.  This  was  detected  by  the  extractor 
and  was  interpreted  as  two  sets  of  comb  actuators.  The  layout 
was  then  corrected  to  remove  the  error.  This  example  demon¬ 
strates  how  the  extractor  can  be  combined  with  visual  inspec¬ 
tion  to  act  as  a  LVS  tool  for  MEMS  designs. 

C.  Gyroscope 

The  issues  like  symmetry  and  cross  axis  coupling  become 
extremely  important  in  complex  devices  like  the  three-fold 
symmetric  gyroscope  [36]  shown  in  Fig.  11(a).  The  gyroscope 
shown  uses  U-springs  and  beams  for  its  suspension  mechanism 
and  pedestal  type  fingers  in  its  comb  drive  for  increased  actu¬ 
ation.  The  extractor  is  able  to  correctly  recognize  the  various 
mechanical  and  electromechanical  parts  of  the  device  resulting 
in  an  extracted  schematic  representation  [see  Fig.  11(b)]  which 


D.  Orthogonal  Accelerometers 

Fig.  12(a)  shows  a  layout  consisting  of  two  orthogonal  ac¬ 
celerometers.  Such  a  combination  is  sometimes  used  to  sense 
acceleration  along  the  orthogonal  axes  [19].  Fig.  12(b)  shows 
the  recognized  representation.  As  can  be  seen,  the  optimization 
algorithm  for  generating  minimal  number  of  plates  and  anchors 
in  the  functional  level  representation  selected  the  correct  rep¬ 
resentation  for  each.  Thus  for  the  accelerometer  on  the  left  it 
selected  the  maximal  horizontal  representation  while  for  the  ac¬ 
celerometer  on  the  right  it  selected  a  maximal  vertical  represen¬ 
tation. 

E.  Experimental  Verification 

An  array  of  folded  flexure  resonators  [35]  were  fabricated  in 
MUMPS  [24]  and  the  experimental  results  [37]  were  compared 
with  the  results  from  simulation  of  extracted  schematic  (using 
NODAS  [6]).  Fig.  13(a)  shows  an  SEM  picture  of  a  30  kHz 
resonator  fabricated  in  the  array.  The  corresponding  extracted 
schematic  is  shown  in  Fig.  13(b).  The  extracted  schematic  in¬ 
cluded  a  mean  overetch  value  of  0.135  pm  which  was  obtained 
from  actual  measurements  of  the  fabricated  structures.  For 
this  particular  example  the  spring  recognition  was  disabled, 
resulting  in  a  schematic  which  contained  atomic  as  well  as 
functional  elements.  This  was  done  to  capture  the  effect  of 
overetch  in  the  beams  which  was  not  available  in  the  functional 
level  model  of  the  folded  flexure  spring  [31],  Accuracy  of 
simulation  result  was  considered  more  important  than  the  gain 
in  simulation  time  that  could  have  been  achieved  if  a  schematic 
with  only  functional  elements  was  used.  The  table  in  Fig.  13(c) 
compares  the  resonant  frequency  and  Q-factor  obtained  from 
simulation  and  experimental  measurements.  The  deviation  in 


10 


JOURNAL  OF  MICROELECTROMECHANICAL  SYSTEMS,  VOL.  11,  NO.  1,  FEBRUARY  2002 


parameter 

measured 

extracted 

extracted  over 
measured 

resonant 

frequency 

22.7  kHz 

22.9  kHz 

0.9% 

quality 

factor 

25.6 

27.1 

5.9% 

(C) 

Fig.  13.  A  30-kHz  folded  flexure  resonator:  (a)  SEM  picture,  (b)  extracted 
schematic,  and  (c)  comparison  of  results. 


on  the  ordinate  axis  obtained  by  projecting  the  edges  in  h  on  the 
ordinate  (n?()  will  also  be  O(r).  If  a  grid  is  formed  using  these 
points  (nv  and  nu)  then  total  number  of  grid  rectangles  ( ng ) 
will  be  0(nv  *  riff)  which  is  Of'2).  Since  canonical  representa¬ 
tion  is  obtained  by  extending  the  outer  edges  of  the  polygon,  the 
number  of  rectangles  in  the  canonical  representation  ( n )  will  be 
less  than  or  equal  to  ng.  Hence  n  has  an  upper  bound  of  Otr2). 

Proposition  2:  The  running  time  for  the  canonization  algo¬ 
rithm  is  Ofnf'a)  where  n  is  the  number  of  rectangles  in  the  final 
canonical  representation. 

Proof:  Total  time  in  canonization  is  the  sum  of  the  time 
required  to  create  the  rectangles  of  the  canonical  representation 
and  the  time  required  to  find  the  ADJACENT  set  (())  in  each 
loop.  The  time  required  to  create  n  rectangles  is  Off).  To  find 
out  the  time  required  to  find  the  ADJACENT  set  let  us  consider 
the  zth  rectangle  of  the  starting  set  11.  The  number  of  rectangles 
already  in  the  final  set  G  at  this  stage  will  be  of  0((i  —  l)2)  [see 
Proposition  1].  Hence,  number  of  comparisons  needed  to  get  Q 
for  the  zth  rectangle  will  be  ()((i  —  l)2).  Summing  up  for  all  i 
from  1  to  r  (where  r  is  the  number  of  rectangles  in  the  starting 
set  R),  we  get,  total  time  in  comparisons  to  be  0(£(z  —  l)2) 
which  is  Of'3)  or  Ofi1^)  [since  n  —  Of' 2)  from  Proposition 
1],  Hence  total  time  in  canonization  is  Ofn  +  n1-°)  which  is 
Ofn1-0). 


Q-factor  is  due  to  the  inaccuracies  in  the  damping  models  used 
in  simulation. 

VI.  Conclusion 

Verification  of  complex  MEMS  designs  need  capabilities  to 
reconstruct  circuit  schematic  representation  from  the  layout  rep¬ 
resentation  of  the  component.  Such  reconstructed  schematics 
can  be  used  to  identify  design  problems  without  performing 
time  consuming  numerical  simulation.  An  extractor  based  on 
geometrical  heuristics  has  been  proposed  and  its  capabilities 
demonstrated.  The  extractor  also  gives  the  user  the  flexibility 
to  extract  to  two  different  levels  of  design  hierarchy  depending 
on  the  user’s  requirements  of  simulation  time  and  accuracy.  The 
prototype  implementation  of  the  extractor  is  limited  to  Man¬ 
hattan  designs  using  the  MUMPS  process,  but  can  be  extended 
to  handle  non  Manhattan  designs  and  other  processes. 

Appendix 

Proposition  1:  The  maximum  size  of  the  number  of  rectan¬ 
gles  (n)  in  the  canonical  representation  of  a  layout  is  Otr2) 
where  r  is  the  number  of  nontrivial  rectangles  used  to  repre¬ 
sent  the  layout. 

Proof:  Let  r  be  the  initial  number  of  nontrivial  rectan¬ 
gles.  Since  the  rectangles  are  nontrivial,  there  are  no  overlap¬ 
ping  rectangles  and  also  no  two  neighboring  rectangles  have  to¬ 
tally  overlapping  edges.  Then,  number  of  vertical  outer  edges  of 
the  layout  (v)  is  Of')  and  similarly,  number  of  horizontal  outer 
edges  of  the  layout  (/?.)  is  also  Of').  Let  nv  be  the  points  on 
the  abscissa  obtained  by  projecting  the  vertical  outer  edges  on 
the  abscissa.  Since  some  of  the  vertical  outer  edges  will  be  co- 
linear,  nv  <  v.  Hence,  nv  is  Of').  Similarly,  number  of  points 
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